ON THE SINGULAR SPECTRUM FOR ADIABATIC QUASI-PERIODIC 
SCHRODINGER OPERATORS ON THE REAL LINE 



ALEXANDER FEDOTOV AND FREDERIC KLOPP 



Abstract. In this paper, we study spectral properties of a family of quasi-periodic Schrodinger op- 
erators on the real line in the adiabatic limit. We assume that the adiabatic iso-energetic curves are 
extended along the momentum direction. In the energy intervals where this happens, we obtain an 
^ - asymptotic formula for the Lyapunov exponent, and show that the spectrum is purely singular. 
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' Resume. Get article est consacre a I'etude du spectre d'une certaine famille d'equations de Schrodinger 

^S) , quasi-periodiques sur I'axe reel lorsque les courbes iso-energetiques adiabatiques sont non bornees 

dans la direction des moments. Dans des intevralles d'energies oil cette propriete est verifiee, nous 
O l' obtenons une formule asymptotique pour I'exposant de Lyapunov, et nous demontrons que le spectre 

, est purement singulier. 

in 

<N . 

0. Introduction 

> 

1/-^ ' In this paper, we continue our analysis of the spectrum of the ergodic family of Schrodinger equations 

cn ■ 

O. ,2 

(0.1) H,^,i; = -^^|J{x) + {V{x-z) + W{ex))i;{x)=E^lJ{x), xe 

m ■ 
o 



where V{x) and VF(^) are periodic and real valued, z € M indexes the equations, and e > is chosen 
so that the potential V{- — z) + W{e-) be quasi-periodic. We study the spectral properties of the 
Q_^| operator H^^e acting in L^(M) in the limit as e — > 0. In the paper [8], we studied this operator near 
the bottom of the spectrum when W is the cosine. In the paper [9], for a general analytic, periodic 
potential W, we studied the spectrum located in the "middle" of a spectral band of the "unperturbed" 
periodic operator 

(0.2) Ho^Pix) = -^"{x) + V{x)^p(x). 
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In the present paper, we again consider a rather general analytic potential W; we only assume that it 
has exactly one maximum and one minimum in a period, and that these are non-degenerate. As about 



■ V, it can be rather singular; for the sake of simplicity, we assume that it belongs to Lf^^. We study 
the spectrum in an energy interval J such that, for all E £ J, the interval E — W{M) contains one or 
more isolated spectral bands of the periodic operator (0.2) whereas the ends of the interval E — W(M) 
are in the gaps, see Fig. 1. So, we are interested in the spectrum close to and inside relatively small 
bands of the unperturbed periodic operator Hq. 

As in [8, 9], our main tool is the monodromy matrix. Most of the present paper is devoted to the 
asymptotic study of the monodromy matrix for the family of equations (0.1). In the adiabatic limit 
e ^ 0, the monodromy matrix is asymptotic to a trigonometric polynomial; if the interval E — W(M) 
contains only one isolated spectral band, this is a trigonometric polynomial of a first order. In result, 
the analysis of (0.1) reduces to the analysis of a "simple" model difference equation. 

Using the monodromy matrix asymptotics, we obtain asymptotic formulae for the Lyapunov exponent 
for the equation family (0.1). They show that, in J, the energy region we study, the Lyapunov 
exponent is positive. This implies that the spectrum of (0.1) in J is singular. 
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Figure 1: The isolated band: two possible cases 



The spectral results admit a natural semi-classical interpretation. Let £-{k) be the dispersion relation 
associated to Hq. Consider the real and the complex iso- energy curves Fk and T defined by 

(0.3) Tr: £{k) + W{Q) = E, k,CgM, 

(0.4) T : S{k) + W{Q) = E, k,CgC. 

These curves are 27r-periodic as in C, so in k. Under our assumptions, the real branches of T (the 
connected components of Tr) are isolated continuous curves periodic in k. In the case when the 
interval E — Ty(M) contains only one spectral band, the iso-energy curve is shown in Fig. 2. The real 
branches are represented by full lines. They are connected by complex loops (closed curves) lying on 
F; the loops are represented by dashed lines. 

The adiabatic limit can be regarded as a semi-classical 
limit, and the expression £{k) + W{C,) can be interpreted as 
a "classical" Hamiltonian corresponding to the operator (0.1). 
Then, from the quantum physicist point of view (see [18, 19]), 
a semi-classical particle should "live" near the real branches of 
the iso-energy curve. In our case, these curves are "extended" in 
momentum variable and "localized" in position variable. There- 
fore, they have to correspond to localized states. The decay 
of these states in the position variable is characterized by the 
complex tunneling between the real branches along the complex 
loops. So, the Lyapunov exponent is naturally related to the 
tunneling coefficients. Our results justify this heuristics. 

This naturally leads to the following conjecture: in a given en- 
ergy interval, if the iso-energy curve has a real branch that is an Figure 2: The phase space picture 
unbounded vertical curve, then, in the adiabatic limit, in this 
interval, the Lyapunov exponent is positive and the spectrum is singular. 

Note that, in [9], we have proved a dual result for the absolutely continuous spectrum: we have proved 
that, if, in some energy region, the branches of the real iso-energy curve are unbounded horizontal 
curves, then, this energy region, except for a set of exponentially small measure, is in the absolutely 
continuous spectrum. 

1. The results 
We now state our assumptions and results. 

1.1. Assumptions on the potential. About the functions V and W , we assume that 
(H): • V and W are periodic, 

(1.1) V {x + l) = V {x), W {x + 2ti) = W {x), xeM; 

• V \s real valued and locally square integrable; 

• W \s real analytic in a neighborhood of M, say, in the strip {|Im2;| < y}; 

• W has exactly one maximum and one minimum in [0, 2tt); they are non degenerate. 




To fix notations, assume that, on the interval [0, 2tt), W is maximum at and minimum at C*- 

In (0.1), e is a positive parameter. For each fixed e, we consider (0.1) as a family of equations indexed 
by the parameter z G M. 

Note that, if 27r/e Q, the function V{x — z) + W{ex) is quasi-periodic in x, the ratio of the frequencies 
of V and W being equal to 27r/e; hence, (0.1) is an ergodic family of equations (see [15]). 

1.2. The assumption on the energy region. To describe the energy regions where we study the 
spectral properties of the family of equations (0.1), we consider the periodic Schrodinger operator Hq 
acting in ^^(M) defined by (0.2). 

1.2.1. Periodic operator. The spectrum of (0.2) is absolutely continuous and consists of intervals of 
the real axis [E'l, E2], [E^, E4], . . . , [S2n+i, -E'2n+2], ■ ■ ■ , such that 

El < E2 < E3 < E4 . . . E2n < E2n+1 < E2n+2 < ■ • • , 

En — > +00, n — > +00. 

The points Ej, j = 1,2,3 ... , are the eigenvalues of the differential operator (0.2) acting on L^{[0, 2]) 
with periodic boundary conditions. The intervals defined above are called the spectral bands, and the 
intervals {E2, E3), [E^, E^), . . . , {E2n, E2n+i), • • ■ , are called the spectral gaps. If E2n < E2n+i, we 
say that the nth gap is open, and, if [£'2^-1, -^2?!] is separated from the rest of the spectrum by open 
gaps, we say that the n-th band is isolated. 

1.2.2. The "geometric" assumption. Let us now describe the energy region where we study the family 
of equations (0.1). 

The spectral window centered at E is the interval W{E) = E — W{M.). If W+ = maXx£M.W{x) and 
W- = min^6M W{x), then, W{E) = [E - W+, E - W.]. 

We assume that there exists J C M, a compact interval such that, for all E ^ J , the window 'W{E) 
contains exactly m + 1 isolated bands of the periodic operator. That is, we fix two integers n > and 
m and assume that 

(Al): the bands [£'2(n+j)-i) -E'2(„+j))], j = 0, 1, . . . m, are isolated; 

(A2): for all E £ J, these bands are contained in the interior of W{E); 

(A3): for all E £ J, the rest of the spectrum of the periodic operator is outside }V{E). 

Note that energies E satisfying (Al) - (A3) exist only if W+ — W-, the "amplitude" of the adiabatic 
perturbation, is large enough; e.g., if m = 0, such energies exist if and only if W+ — W- is larger than 
the size of the n-th spectral band, but smaller than the distance between the (n — l)-st and (n -|- l)-st 
bands. 

From now on, unless stated otherwise, we assume that our assumptions on V and W, and assumptions 
(Al) - (A3) are satisfied. 

1.3. Iso-energy curve. Our results are formulated in terms of the iso-energy curve F defined by (0.4). 
The iso-energy curve is 27r periodic both in the ^— and k— directions (see Lemma 10.1). 

1.3.1. The real branches. To describe the real branches of F, i.e. the connected components of the real 
iso-energy curve Fjr, we define the following collection of subintervals of [0, 2tt]. Consider the mapping 

£:C^E-W{0. 

It is monotonous on each of the intervals /_ = [0, ("*] and /+ = [C*, 2tt] and maps each of them onto 
the spectral window W{E). For n < j < n + m, let 3^ C /+ (resp. C /_) be the the pre-image 
of the j-ih. spectral band in W{E). Let Z be the collection of these intervals. A "period" of the real 
iso-energy curve is described by 

Lemma 1.1. Let E £ J. The set F^ n {0 < k < 27r} consists of 2{m + 1) curves {7(3), 3 G Z}. 
Fix I £ Z. The curve 7(3) is the graph {(k, (") : C = ^lii^)-: At G 1^} of a function which satisfies 

(1) it is continuous, 

(2) it is 2'K-periodic and even in k, 

(3) it is monotonous on the interval [0,7r], 
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(4) it maps [0, vr] onto 3. 
The curves 7(3) continuously depend on E £ J. 

Lemma 1.1 is proved in section 10.1. For m = 0, the real iso-energy curve is shown in Fig. 2. 

1.3.2. Complex loops. Now, we discuss loops, i.e. closed curves, situated on the iso-energy curve P 
and connecting its real branches. 

For J = n — 1, n, . . . , n + m, let 0^ (resp. qJ) be the subinterval of (resp. /_) that is the pre-image 
of the part of j-th spectral gap situated inside W{E). Let 

0n-l = (0^_i - 27r) U 0~_i and Qn+m = et+rn U 0n+m- 

Then, 0n-i is an open interval containing zero, and Qn+m is an open interval containing Let Q be 
the set consisting of Qn-i, Qn+m and the intervals with j = n,n + 1, . . . ,n + m — 1. 

For Q G, let V{q) C C be a sufficiently small complex neighborhood of the interval 0. Let G{q) be 
a smooth closed curve that goes once around the interval in V{q) \ q. In Figure 3, we depicted the 
curves G(0) when m = 0. 

In section 10.2, we show that each of the curves G{q) is the projection of ^^(0), a closed curve on F. 
This curve connects the real branches projecting onto the intervals adjacent to 0. 

1.3.3. Tunneling coefficients. To F, we associate the tunneling coefficients 

(1.2) t(0) = 6-^-^(0), qgG, 
where S{g) are the tunneling actions given by 

(1.3) 5(0) = i I KdC, qgG. 

In section 10.3, we show that, for E £ J, each of these actions is real and non-zero. By definition, we 
choose the direction of the integration so that all the tunneling actions be positive. 
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Figure 3: The curves G(0) for m = 



1.4. Spectral results. One of the main objects of the spectral theory of quasi-periodic equations is 
the Lyapunov exponent, see, for example, [15]. Our main spectral result is 

Theorem 1.1. Let J be an interval satisfying the assumptions (A1)-(A3) for some n and m. Let 
W and V satisfy the hypothesis (H), and let e he irrational. Then, on the interval J , for sufficiently 
small, irrational e/27r, the Lyapunov exponent Q{E) for the family of equations (0.1) is positive and 
has the asymptotics 

(1.4) e(i;)=ij:i„^+o(i)=ij:s(o)+o(i). 

gee gee; 

Note that, this theorem implies that, if e is sufficiently small, then, the Lyapunov exponent is positive 
for all E £ J. 

Recall that, if 2Tr/e is irrational, then Hz,e is quasi-periodic. In this case, its spectrum does not depend 
on z (see [1]); denote it by a{Hz^e)- In l^]-, we have proved 
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Theorem 1.2 ([S]). Let S = a{Ho) + W{R) = a{Ho) + W+]. Then, one has 

• Ve > 0, cr{H,^e) C S. 

• for any K C T, compact, there exists C > such that for all e sufficiently small and \/E G K , 
one has 

a{H,,e) r\{E- E + Ce^/^) / 0. 

By the Ishii-Pastur-Kotani Theorem [4, 15] and Theorem 1.5 in [ Theorems 1.1 and 1.2 imply 
Corollary 1.1. In the case of Theorem 1.1, for e sufficiently small, for all z G M, one has 

a{H^^e) nJ^$ and aac{H^^e) r\J = ^, 

where o'ac{Hz,e) the absolutely continuous spectrum of the family of equations (0.1). 

1.5. The monodromy matrix and Lyapunov exponents. The main object of our study is the 
monodromy matrix for the family of equations (0.1); we define it briefly (we refer to [7, 8] for more 
details). The central result of the paper is its asymptotics in the adiabatic limit. 

1.5.1. Definition of the monodromy matrix. Consider a consistent basis (V'1,2) i-e. a basis of solutions 
of (0.1) whose Wronskian is independent of z and that are 1-periodic in z i.e. that satisfy 

(1.5) ipi^2{x, z + I) = ipi^2ix, z), Mx,z. 

The functions ^'1,2(2^ + 27r/e, z + 27r/e) being solutions of equation (0.1), one can write 

(1.6) ^{x + 27r/e, z + 27r/e) = M {z) ^ {x, z), 
where 

. ^(x z) = f'^i^^' 

• M (z) is a 2 X 2 matrix with coefficients independent of x. 

The matrix M is called the monodromy matrix associated to the consistent basis (V'1,2 )• Note that 

(1.7) detM{z) = l, M{z + 1) = M{z), \/z. 

1.5.2. Monodromy equation and Lyapunov exponents. Set h = ^modl. Let M be the monodromy 
matrix associated to a consistent basis (V'1,2)- Consider the monodromy equation 

(1.8) Fn+i = M{z + nh)Fn Vn G Z. 

There are several deep relations between the monodromy equation and the family of equations (0.1) 
(see [9, 8]). We describe only one of them. Let 27r/e be irrational, and let Q{E) (resp. 0{E)) be the 
Lyapunov exponent for (0.1) (resp. for (1.8)). One proves 

Theorem 1.3 {[>]). The Lyapunov exponents Q(E) and 9{E) satisfy the relation 

(1.9) Q{E) = ^e{E). 

The passage to the monodromy equation is close to the monodromization idea developed in [2] for 
difference equations with periodic coefficients. 

1.5.3. The asymptotics of the monodromy matrix. As W and V are real on the real line, we construct 
a monodromy matrix of the form 



(1.10) 



a{z,E) b{z,E) 
b{z,E) a{z,E) 



In the adiabatic case, the asymptotics of a and b have very simple, model form. We first assume that 
n in (Al) - (A3) is odd. Then, one has 
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Theorem 1.4. Let E € J. There exists y > and Vq, a neighborhood of Eq, such that, for sufficiently 
small e, the family of equations (0.1) has a consistent basis of solutions for which the corresponding 
monodromy matrix M is analytic in {z,E) G {|Imz| < Y/e] x Vq and has the form (1.10). The 
coefficients a and b admit the asymptotic representations 

(1.11) a = a_^e-2'^^™^(l + o(l)), 6 = e-2^^"*^(l + o(l)), 0<lmz<Y/e, 
and 

(1.12) a = a„+ie2-('"+^)-(l + o(l)), b = bm+i e^^'^^+^> {1 + o{l)) , -Y/e<lmz<0. 

The coefficients a^jji, fe— mi o,m+i and bm+i fffi independent of z. Moreover, there exists a constant 
C > 1 ( independent of e and E ) such that 

^ <T ■\aj\<C, ^ <T ■\bj\<C, where T = T{E) = \[^^g t(g), 

(1.13) 

G Vb n M, J = -m, m + l. 

Pick Yi and Y2 so that < Yi < I2 < There is V = V{Yi,Y2), a neighborhood of Eq such that the 
asymptotics of a and b are uniform in {z,E) £ {Yi < |Im(^| < I2} x ^■ 

In sections 9.1 and 9.2, we give asymptotic formulae for Om+i, o--m and bm+i, b^rn- 

In the case n even, one has a similar result. The only novelty is that, in this case, the formulae (1.11) 
and (1.12) describe the asymptotics of the coefficients of the matrix related to M, the monodromy 
matrix, by the following transformation 

(1.14) S-^{z + h)M[E,z)S{z), S{z) = 

The asymptotics (1.11) and (1.12) are obtained by means of the new asymptotic method developed 
in [7, 10]. 

1.5.4. Fourier coefficients. The coefficients Om+i, o-m and bm+i, fe-m are the leading terms of the 
asymptotics of the (m + l)-th and (— m)-th Fourier coefficients of the monodromy matrix coefficients. 
Theorem 1.4 implies that, in the strip {|Im<^| < y}, the leading terms of the asymptotics of the 
monodromy matrix are equal to the contribution of a few of its Fourier series terms. 

1.5.5. The case m = 0. When m = (and n odd). Theorem 1.4 imply that, in the whole strip 
{|Imz| < Y}, the monodromy matrix coefficients a and b admit the asymptotics: 

(1.15) o = ao (1 + 0(1)) + ai e^^'^'il + o(l)), 6 = 60 (1 + o(l)) + 61 e^'^'^il + o(l)). 

So, up to the error terms, the monodromy matrix becomes a first order trigonometric polynomial: 

(1.16) M~Mo=f^°+^^" bo + hu\ 2... 
^ ^ \bo + bi/u ao + ai/uj' 

with constant coefficients qq, oi, bo, bi of order 0{1/T{E)) (for real E). 

We see that, for m = the monodromy equation becomes a "simple" model equation. 

1.5.6. Relation to the spectral results. In this paper, we use the asymptotics of the monodromy matrix 
only to prove Theorem 1.1. However, we believe that these asymptotics can be used to get quite a 
detailed information on the spectrum of (0.1) in the adiabatic limit. Therefore, we plan to study 
the model equation with the matrix Mq in a subsequent paper. In particular, it seems reasonable 
to believe that, under a Diophantine condition on 27r/e, the spectrum of (0.1) is pure point and the 
eigenvalues can be described by quantization conditions of Bohr-Sommerfeld type. 
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1.5.7. Organization of the paper. Section 2 is devoted to the proof of Theorem 1.1 using Theorem 1.4. 
In section 3, we recall some well known facts from the theory of periodic Schrodinger operators on the 
real line. In section 4, we recall the main construction of the asymptotic method we use to compute 
the monodromy matrix. In sections 5 and 6, using this method, we construct a consistent basis of 
solutions having a simple "standard" asymptotic behavior in the complex plane of ( = ez. In section 7, 
we discuss the properties of the monodromy matrix for this basis. This is the monodromy matrix the 
asymptotics of which are described in Theorem 1.4. Sections 8 and 9 are devoted to the computation 
of the asymptotics of the monodromy matrix. In section 10, we study the geometry of the iso-energy 
curve r and prove estimates (1.13). 

2. The asymptotics for the Lyapunov exponent 

In this section, we prove the asymptotics (1.4). We deduce these asymptotics from the asymptotics 
of the monodromy matrix coefficients described by Theorem 1.4. First, we use a statement of [S] 
and obtain a lower bound for the Lyapunov exponent. This statement is based on the ideas of [IG] 
generalizing Herman's argument [1 1] . Then, using the asymptotics of the coefficients of the monodromy 
matrix in the complex plane, we get estimates on the real line. This yields an upper bound for the 
Lyapunov exponent. Comparing the upper and the lower bounds, we obtain (1.4). 

Recall that, for equation (1.8), the Lyapunov exponent is defined by 

(2.1) 9{M) = lim llog||Pjv(^)||, 

where Pn is the matrix cocycle 

Pn{z) = M{z + Nh) ■ M{z + {N- l)h) • • • M{z + h) ■ M{z). 

It is well known (see [3, 15, 16] and references therein) that, if h is irrational, and M{z) sufficiently 
regular in z, then the limit (2.1) exists for almost all z and is independent of z. 

2.1. The lower bound. 

2.1.1. Preliminaries. Let (M(2;, e))o<e<i be a family of S'L(2, C)-valued 1-periodic functions of z G C. 
Let h be an irrational number. One has 

Proposition 2.1 ([S]). Pick eq > 0. Assume that there exist yo ond yi satisfying the inequalities 
< t/o < yi < oo and such that, for any e G (0, eo) one has 

• the function z M{z,£) is analytic in the strip S = {z G C; < Imz < yi/e}; 

• in the strip Si = {z ^ C; ya/e < Imz < yi/e} C S, M{z,£) admits the representation 

M{z,e) = A(e)e^2™«^ • {Mo{e) + Mi{z,e)) , 
for some constant A(e), an integer no and a matrix Mo(e), all of them independent of z; 

• there exist constants /3 > and a G (0,1) independent of e and such that |a(e)| < a and 

\m\<P; 

• sup^gs^ ||Mi(2;,e)|| < m{e), m{e) ^0 as e ^ 0. 

Then, there exit C > and ei > (both depending only on y^, yi, a, (3 and m{-)) such that, if 
< e < ei, one has 

(2.2) e{M) > log |A(e)| - Cm(e). 

In [8], we have assumed that no is a positive integer, but the proof remains the same for no G Z. 

We use this result and Theorem 1.4 to get the lower bound for the Lyapunov exponent. In the sequel, 
n is the index introduced in assumptions (Al) - (A3). The cases n odd and n even are treated 
separately. 
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2.1.2. Obtaining the lower bound for n odd. In the sequel, we suppose that the assumptions of The- 
orem 1.4 are satisfied; we use the notations and resuhs from this theorem without referring to it 
anymore. We assume that E e Vq HM.. 

'0 r 



Let <^ = ^-|^ gj ■ Show that the matrix aM{z)a satisfies the assumptions of Proposition 2.1. 

Fix 2/0 and yi so that < jjq < yi < Y . The asymptotics of the monodromy matrix coefficients are 
uniform in z in the strip S = {yo/e < Imz < yi/s} and hi E £ Vq (reducing Vq if necessary). For 
E £ Vq CiM and z £ S, formulae (1.11) and (1.12), and estimates (1.13) imply that 



aiz) = a;;^e-2-(-+i)^(l + o(l)), M = e{E) (1 + o(l)), ^ = o(l), =^ = o(l) 
where c{E) is independent of z and bounded by a constant uniformly in e and E. So, we have 



a M{z)a = Um+i e 



-2TTi(m+l)z 



1 c{E) 




+ o(l] 



We see that the matrix aM{z)a satisfies the assumptions of Proposition 2.1. 

Clearly, the Lyapunov exponents of the matrix cocycles associated to the pairs (M, h) and {aMa, h) 
coincide. So, Proposition 2.1 implies that 0{M), the Lyapunov exponents of the matrix cocycle 
associated to (M, h), satisfies the estimate 0{M) > log |a_m| + o(l)- 

The Lyapunov exponent Q{E) for equation (0.1) is related to 6{M) by Theorem 1.3. Therefore, 
&iE) > ^log|a_m| +o{e). Hence, (1.13) clearly implies 

(2.3) e(i?)>^iogr-i + o(e). 

2.1.3. The lower bound when n is even. If n is even, then, formulae (1.11) and (1.12) give the asymp- 
totics of the coefficients of the matrix (1.14). Obviously, the Lyapunov exponents for the matrix 
cocycles generated by M{z) and by S~^{z + h)M{z)S{z) coincide. Arguing exactly as in subsec- 
tion 2.1.2, we again obtain (2.3). 

2.2. The upper bound. Let us first assume that n is odd. Let E £ Vq DM.. Fix < yo < The 
asymptotics (1.11) and (1.12) and estimates (1-13) imply the following estimates for the coefficients 
of M{z), the monodromy matrix: 

|a|,|6| < C(yo)T-ie2™W^, Imz = yo/e, 
|a|,|6| < C(?/o)^-le2'^('"+l)^^o/^ Imz = -yo/e. 

Here, C{yo) is a positive constant independent of e, Re z, and E. The estimates are valid for sufficiently 
small e. Recall that M is analytic and 1-periodic in z. Therefore, (2.4) and the Maximum Principle 
imply that 

(2.5) \a\,\b\ <2C{yo)T-'^exp{2TT{m + l)yo/e), z£R. 

This leads to the following upper bound for the Lyapunov exponent for the matrix cocycle generated 
by M 

9 < logT"^ + Const + 27r(m + l)yo/e 
where Const is independent of E and e. In view of Theorem 1.3, we finally get 

(2.6) @{E)<— log + eConst + 27r(m + l)yo- 

The upper bound (2.6) remains true when n is even as the Lyapunov exponents for the matrix cocycles 
generated by M{z) and by S^^{z + h)M{z)S{z) coincide. 

2.3. Completing the proof. Recall that, in (2.6), yo is an arbitrarily fixed positive number. So, 

comparing (2.3) and (2.6), we get Q{E) = — logT"^ -|- o(l). This and the formula (1.13) for T 

27r 

imply (1.4) for ah E£VonR. 

Recall that Vq H M is an open interval containing Eq £ J. The above construction can be carried out 
for any Eo £ J. As the interval J is compact, this completes the proof of Theorem 1.1. 



3. Periodic Schrodinger operators 



We now discuss the periodic Schrodinger operator (0.2) where F is a 1-periodic, real valued, Lf^^- 
function. We collect known results needed in the present paper (see [5, 13, 14, 17, 10]). 

3.1. Bloch solutions. Let be a solution of the equation 

(f 

(3.1) -—.i){x) + Vix)'4){x)=8^{x), xGM, 

satisfying the relation ij: [x + l) = \ip {x) for all x G M with A G C independent of x. Such a solution is 
called a Bloch solution^ and the number A is called the Floquet multiplier. Let us discuss the analytic 
properties of f i— > ^p{S, x) := ip{x). 

As in section 1.2, we denote the spectral bands of the periodic Schrodinger equation by [Ei, E2], 
[E3, E4], . . . , [E2n+i, E2n+2], ■ ■ ■ ■ Consider S±, two copies of the complex plane £ £ C cut along the 
spectral bands. Paste them together to get a Riemann surface with square root branch points. We 
denote this Riemann surface by S. 

One can construct a Bloch solution ip{x,£) of equation (3.1) meromorphic on S. It is normalized by 
the condition ^{1,£) = 1. The poles of this solution are located in the spectral gaps. More precisely, 
each spectral gap contains precisely one simple pole. It is located either on 5+ or on The position 
of the pole is independent of x. 

For £ G S, we denote by £ the point on S having the same projection on C as £. We let 

iIj{x.,£) = il){x,£), ££S. 

The function ij){x,£) is another Bloch solution of (3.1). Except at the edges of the spectrum (i.e. the 
branch points of 5), the functions ij: and are linearly independent solutions of (3.1). In the spectral 
gaps, "0 and il) are real valued functions of x, and, on the spectral bands, they differ only by complex 
conjugation. 

3.2. The Bloch quasi-momentum. Consider the Bloch solution Tp{x,£). The corresponding Flo- 
quet multiplier A {£) is analytic on S. Represent it in the form X{£) = exp{ik{£)). The function k(£) 
is the Bloch quasi-momentum. 

The Bloch quasi-momentum is an analytic multi-valued function of £. It has the same branch points 
as tpix, £). 

Let Z) be a simply connected domain containing no branch point of the Bloch quasi-momentum. In 
D, one can fix an analytic single- valued branch of k, say kQ. All the other single- valued branches of k 
that are analytic in D are related to /cq by the formulae 

(3.2) k±^iiS) = ±koi£) + 2ttI, I G Z. 

Consider C+ the upper half plane of the complex plane. On C-|-, one can fix a single valued analytic 
branch of the quasi-momentum continuous up to the real line. It can be fixed uniquely by the condition 
—ik{£ + iQ) > Q as £ < El. We call this branch the main branch of the Bloch quasi-momentum and 
denote it by kp. 

The function kp conformally maps C+ onto the first quadrant of the complex plane cut at compact 
vertical slits starting at the points vr/, / G N. It is monotonically increasing along the spectral zones 
so that [£'2„_i, -E'2n]) tlie n-th spectral band, is mapped on the interval [7r(n — l),7rn]. Along any 
open gap, Ylekp{£ + iO) is constant, and Imfep(£^ -|- iO) is positive and has only one non-degenerate 
maximum. 

All the branch point of kp are of square root type. Let Ei be a branch point. In a sufficiently small 
neighborhood of Ei, the function kp is analytic in yj£ — Ei, and 

(3.3) kpi£) - kp{Ei) = ciy/£-Ei + 0{£ - Ei), ci^O. 

Finally, we note that the main branch can be analytically continued on the complex plane cut only 
along the spectral gaps of the periodic operator. 
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3.3. A meromorphic function. Here, we discuss a function playing an important role in the adia- 
batic constructions. 

In [10], we have seen that, on S, there is a meromorphic function uj having the following properties: 

• the differential = lu d£ is meromorphic; its poles are the points of P U Q, where P is the set 
of poles of Tp{x,£)^ and Q is the set of points where k'{£.) = 0. 

• all the poles of are simple; 

• mSpQ. = l, yp£P\Q, reSgJ^ = -l/2, G Q \ P, res ^f) = 1/2, VrGPnQ. 

• if £ & S projects into a gap, then co{£) G M; 

• ii £ £ S projects inside a band, then lo{£) = uj{£). 



4. The complex WKB method for adiabatic problems 

In this section, following [7, 10], we briefly describe the complex WKB method for adiabatically 
perturbed periodic Schrodinger equations 

(4.1) --^Tp{x) + {V{x) + W{ex + C))il:{x) = EiPix), x£R. 

Here, V is 1-periodic and e is a small positive parameter; one assumes that V G Lf^^ and that W is 
analytic in T>{W), a neighborhood of the real line {W is not necessarily periodic). 

The parameter is an auxiliary parameter used to decouple the "slow variable" = ex and the "fast 
variable" x. The idea of the method is to study solutions of (4.1) on the complex plane of C and, 
then to recover information on their behavior in x G M. Therefore, one studies solutions satisfying the 
condition: 

(4.2) ^^a: + l,C) = Hx,C + e) VC 

On the complex plane of C, there are certain domains on which one can construct such solutions that, 
moreover, have simple asymptotic behavior. 

4.1. Standard behavior of solutions. We first define two analytic objects central to the complex 
WKB method, the complex momentum and the canonical Block solutions. Then, we describe the 
standard behavior of the solutions studied in the framework of the complex WKB method. 

4.1.1. The complex momentum. The complex momentum k is the main analytic object of the complex 
WKB method. For ^ G P(VF), the domain of analyticity of the function W, it is defined by the formula 



(4.3) K{C) = k{£{0), £{0 = E-W{C), 

Here, k is the Bloch quasi-momentum of (0.2). Relation (4.3) "translates" properties of k into prop- 
erties of K. Hence, k is a multi-valued analytic function, and that its branch points are related to the 
branch points of the quasi-momentum by the relations 

(4.4) Ei + W{C) = E, / = 1,2,3... 

Let ^0 be a branch point of k. If W'{Co) Oi then is a branch point of square root type. 

If D C ^{W) is a simply connected set containing no branch points of k, we call it regular. Let Kp be 
a branch of the complex momentum analytic in a regular domain D. All the other branches analytic 
in D are described by the formulae: 

(4.5) = ±Kp + 27rm, 
where it and ?n G Z are indexing the branches. 
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4.1.2. Canonical Block solutions. To describe the asymptotic formulae of the complex WKB method, 
one needs Bloch solutions to the equation 

(4.6) -^^{x) + V{x)i:{x)=£{CMx), £{C) = E-W{0, xeR, 
that, moreover, are analytic in on a given regular domain. 

Pick (^0 a regular point. Let £q = £{Co)- Assume that £q ^ P U Q. Let Uq be a sufficiently small 
neighborhood of Sq, and let Vq be a neighborhood of such that £{Vo) C Uq. In Uq, we fix a branch 
of the function \jk!{£) and consider '\\)±{x,E\ two branches of the Bloch solution ip(x,S), and 
the corresponding branches of fi. Put 

(4.7) ^i{x,C)=q{£)e^^o''^^^J±ix,£), q{£) = ,/¥{8), £ = £{C). 

The functions ^± are called the canonical Bloch solutions normalized at the point (^q. 

The properties of the differential Q imply that the solutions ^± can be analytically continued from 
Vq to any regular domain D containing Vq. 
One has (see [10]) 

(4.8) w{^+{; C), ^-i; 0) = w{^+{; Co), ^-(^ Co)) = k' {£q)w{4;+{x , ^Tq), V'-(x, ^:o)) 

As <So <5 U {El}, the Wronskian u;(^'+(-, C), ^-(-j C)) is non-zero. 

4.2. Solutions having standard asymptotic behavior. Fix E = Eq. Let D be a regular domain. 
Fix Co G so that if (Co) ^ P^Q. Let k be a branch of the complex momentum continuous in D, and 
let ^'-1- be the canonical Bloch solutions defined on D, normalized at Co and indexed so that n be the 
quasi- momentum for ^'_|_. 

Definition 4.1. Let cr be either -|- or — . We say that, in D, a solution / has standard behavior (or 
standard asymptotics) / ~ exp(af /^KdC) • if 

• there exists Vq, a complex neighborhood of Eq, and X > Q such that / is defined and satis- 
fies (4.1) and (4.2) for any {x,C,E) £ [-X,X] x D x Vq; 

• / is analytic in ( £ D and in E £ Vq; 

• for any K, a compact subset of D, there is ^ C Vq, a neighborhood of Eq, such that, for 
{x, C, E) E [—X, X] X K X V, f has the uniform asymptotic 

(4.9) / = e^l/'^'^^^(vI/, + o(l)), as e ^ 0, 

• this asymptotic can be differentiated once in x without loosing its uniformity properties. 

4.3. Canonical domains. Canonical domains are important examples of domains where one can 
construct solutions with standard asymptotic behavior. They are defined using canonical lines. 

4.3.1. Canonical lines. A curve is called vertical if it is connected, piecewise C^, and if it intersects 
the lines {ImC = Const} at non-zero angles 6, < 6 < tt. Vertical curves are naturally parameterized 
by ImC. 

Let 7 be a regular curve. On 7, fix k, a continuous branch of the complex momentum. 

Definition 4.2. The curve 7 is canonical if it is vertical and such that, along 7, 

(1) Im J*' Kcl^ is strictly monotonously increasing with ImC, 

(2) Im /^(k — 7r)dC is strictly monotonously decreasing with ImC- 

Note that canonical lines are stable under small C^-perturbations. 

4.3.2. Canonical domains. Let K he a, regular domain. On K, fix a continuous branch of the complex 
momentum, say k. The domain K is called canonical if it is the union of curves canonical with respect 
to K and connecting two points Ci and C2 located on dK. 

One has 
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Theorem 4.1 ([10]). Let K be a bounded domain canonical with respect to k. For sufficiently small 
positive e, there exists {f±), two solutions of (4.1), having the standard behavior in K: 



For any fixed x G M, the functions f±{x,C,) ^^'"c analytic in Q in the smallest strip {Yi < ImC, < Y2} 
containing K . 

One easily calculates the Wronskian of the solutions f±{x, C,) to get 



By (4.8), for C, in any fixed compact subset of K and e sufficiently small, the solutions f± are linearly 
independent. 

4.4. The strategy of the WKB method. Our strategy to apply the complex WKB method is 
explained in great detail in [iU]; we recall it briefly. First, we find a canonical line. Roughly, we 
build it out of segments of some "elementary" curves described in section 4.6. Then, we find a "local" 
canonical domain K "enclosing" this line, see section 4.5. For this domain, we construct the solutions 
f± using Theorem 4.1. 

Second, we describe the asymptotic behavior of f± outside the domain K. Therefore, we use three 
general principles presented in section 4.7. 

Having investigated the behavior of f± for —X < x < X and a sufficiently large set of Ci we recover 
the behavior of f± on the real line of x by means of condition (4.2). 

Below, we assume that D is a regular domain, and that k is a branch of the complex momentum 
analytic in D. A segment of a curve is a connected, compact subset of that curve. 

4.5. Local canonical domains. Let 7 C -D be a canonical line (with respect to k). Denote its ends 
by Ci and ^2- Let a domain K d D he a, canonical domain corresponding to the triple k, Qi and ^2- If 
-f £ K, then, K is called a canonical domain enclosing 7. One has 

Lemma 4.1 ([n]). One can construct a canonical domain enclosing any given canonical curve. 
Canonical domains whose existence is a consequence of this lemma are called local. 

4.6. Pre-canonical lines. To construct a local canonical domain we need a canonical line. To con- 
struct such a line, we first build a pre-canonical line made of some "elementary" curves. 

Let 7 C -D be a vertical curve. We call 7 pre-canonical if it is a finite union of segments of canonical 
lines and/or lines of Stokes type i.e. the level curves of the harmonic functions C ^ I™ ndC, or 
C i-^- Im /""(ac — 'K)dC,. One has 

Proposition 4.1 ([''~']). Let ^ be a pre-canonical curve. Denote the ends of j by Qa and (b- 

For V C D, a neighborhood of j and Va C D, a neighborhood of (a, there exists a canonical line 7 C F 

connecting the point ^6 to a point in Va. 

When constructing pre-canonical lines, one uses lines of Stokes type. To analyze them, one uses 

Lemma 4.2. The lines of Stokes type of the family Im ndC, = Const are tangent to the vector field 
k{C); those of the family Im /^(k — vr) = Const are tangent to the vector field k{Q — vr. 

4.7. Tools for computing global asymptotics. A set is said to be constant if it is independent of 




(4.10) 



w{f+,f^) 



ii;(^'+,*_) -ho(l). 



e. 
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4.7.1. The Rectangle Lemma: asymptotics of increasing solutions. The Rectangle Lemma roughly says 
that a solution / preserve the standard behavior along a line Im = Const as long as the leading term 
of the standard asymptotics increases. 

Fix r]m < TjM- Define 5 = {C G C : f?m !^ < ?7m}- Let 71 and 72 be two vertical lines such that 
71 n 72 = 0. Assume that both lines intersect the strip S at the lines Im^ = r/m and Im^ = r^jv/, and 
that 71 is situated to the left of 72 . 

Consider the compact R bounded by 71, 72 and the boundaries of S. Let D=R \ (71 U 72). One has 

Lemma 4.3 (The Rectangle Lemma [!)]). Assume that the ''rectangle" R is contained in a regular 
domain. Let f he a solution to (4.1) satisfying (4.2). Then, for sufficiently small e, one has 

- ndC 

1: If ImK < in D, and if f has standard behavior f ~ in a neighborhood 0/71, 

then, it has standard behavior in a constant domain containing the "rectangle" R. 

2: If Iran > in D, and if f has standard behavior f ~ ^+ in a neighborhood 0/72, 

then, it has the standard behavior in a constant domain containing the "rectangle" R. 

Lemma 4.3 was proved in [9] where one can find more details and references. 

4.7.2. The Adjacent Canonical Domain Principle. This principle complements the Rectangle Lemma; 
it allows us to obtain the asymptotics of decreasing solutions. 

Let 7 be a vertical curve. Let S be the minimal strip of the form {Ci < ImC < C2} containing 7. Let 
U C S he a regular domain such that 7 C dU. We say that U is adjacent to 7. One has 

Proposition 4.2 (The Adjacent Canonical Domain Principle [9]). Let ^ be a canonical line. Assume 
that f, a solution to (4.1) satisfying (4.2) has standard behavior in a domain adjacent to 7. Then, f 
has the standard behavior in any bounded canonical domain enclosing 7. 

4.7.3. Adjacent canonical domains. To apply the Adjacent Canonical Domain Principle, one needs 
to describe canonical domains enclosing a given canonical line. It can be quite difficult to find a 
"maximal" canonical domain enclosing a given canonical line. In practice, one uses "simple" canonical 
domains described in 



Lemma 4.4 (The Trapezium Lemma [10]). Let 70 be a canonical line. 

Triangle 
\7o 



Let U be a domain adjacent to 7, a canonical line containing 7n as rr^ 

A , TT T , T / I .trapezium 

an internal segment. Assume that, m U , ImK 7^ 0. Let (resp. a^j 





be the line of Stokes type beginning at the upper (resp. lower) end of 

7o and going in U downward (resp. upward) from this point. \ 1 \ /^^rf 

One has: 

Trapezium case: Let T C U be a regular domain bounded by Ou, Cd, 70 Oi,nd 7, one more 
canonical line not intersecting 79. Then, T is part of a canonical domain enclosing 79. 

Triangle case: Assume that intersects a^- Let T C U be a regular domain bounded by au, 
ad and line 70. Then, T is part of a canonical domain enclosing 70. 

There always exists a canonical line containing 70; moreover, if Im /t 7^ 
in [/, then, the lines 0"^ and described in the Trapezium Lemma 
always exist (see Lemma 5.3 from [10]). 

4.7.4. The Stokes Lemma. The domains where one justifies the standard behavior using the Adjacent 
Canonical Domain Principle are often bounded by Stokes lines (see definition below) beginning at 
branch points of the complex momentum. The Stokes Lemma allows us to justify the standard 
behavior beyond these lines by "going around" the branch points. 

Notations and assumptions. Assume that Co is a branch point of the complex momentum such that 
W^'(Co) / 0. 

Definition 4.3. The Stokes lines starting at Co are the 
curves 7 defined by 



Im [\niC) - ^{Co))dC = 0, CG7- 
-'Co 
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The angles between the Stokes hnes at Co are equal to 2tt/3. We denote them by cri, cj2 and so 
that (Ti is vertical at Co (see Fig. 4). 

Let o"! be a (compact) segment of ai which begins at Co) is vertical and contains only one branch 
point, i.e. the point Co- 
Let y be a neighborhood of cJi. Assume that V is so small that the Stokes lines fJi, cj2 and divide 
it into three sectors. We denote them by Si, S2 and ^3 so that Si be situated between ai and cj2, and 
the sector ^2 be between a2 and (T3 (see Fig. 4). 

The statement. In [10], we have proved 

Lemma 4.5 (Stokes Lemma). Let V be sufficiently small 
Let f be a solution to (4.1) satisfying (4.2) that 

has standard behavior f ~ -^^ '^'^''^-|_ inside the sector 
5*1 U (72 U ^2 of V. Moreover, assume that, in Si near 
ai, one has ImK(C) > if Si is to the left of ai and 
Im k{Q > otherwise. 

Then, f has standard behavior inside V \cri, the leading 
term of the asymptotics being obtained by analytic contin- 
uation from 5*1 U (T2 U S2 into V' . 

Comments and details on this lemma can be found 
j-^Qj Figure 4-' The Stokes lines in a neighbor- 

hood of a branch point 

5. A CONSISTENT SOLUTION 

We now begin the construction of the consistent basis the monodromy matrix of which we compute. 
Recall that V and W satisfy the hypothesis (H), and fix E = Eq £ J, an interval satisfying the 
hypothesis (Al) - (A3), 

In the present section, by means of the complex WKB method, we construct and study / a solution 
of (0.1) satisfying the consistency condition (1.5). 

To use the complex WKB method, we rewrite (0.1) in terms of the variables 
(5.1) X := X — z, and C = 

It then takes the form (4.1). In the new variables, the consistency condition (1.5) becomes (4.2) 

We first describe the complex momentum and Stokes lines. Then, we construct a local canonical 
domain, hence, a consistent solution to (4.1) by Theorem 4.1. Finally, using the continuation tools, 
we describe global asymptotics of this solution. 

5.1. The complex momentum. We begin with the analysis of the mapping 8 : Q ^ E — W{C,). 
5.1.1. The set W~^{m). As S G M, £:-i(M) coincides with W^^{M). 

The set Ty~"^(M) is 27r-periodic. It consists of the real line and of complex branches (curves) symmetric 
with respect to the real line. There are complex branches separated from the real line, and complex 
branches beginning at the real extrema of W . These do not return to the real line. 

Consider an extremum of W on the real line, say Co- By assumption (H), it is non-degenerate. This 
implies that, near Co; the set Vl^~^(M) consists of a real segment, and of a "complex" curve symmetric 
with respect to the real axis, intersecting the real axis at Co only. This curve is orthogonal to the real 
line at Co- 

For y > 0, we let 5y = {—Y < ImC < Y}. We assume that Y is so small that 

• Sy is contained in the domain of analyticity of W; 

• the set Ty~^(R) H Sy consists of the real line and of the complex lines passing through the real 
extrema of W; 

For such Y, and if m = 0, the set P^~^(M) n Sy is shown in Fig. 5. 
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5.1.2. Branch points. The branch points of the complex momentum are related to the branch points 
of the Bloch quasi-momentum by equation (4.4). So, they lie on VF~"'^(M), and form a 27r-periodic set. 

Consider the branch points situated in the interval [0, 27r) of the real line. Recall that, by assumption, 
the function (" i— > VF(C) has two extrema in [0,27r]: a maximum at C = and a minimum at C*, 
< C* < 27r. The mapping £ is monotonous on each of the intervals /_ = [0, C^,] and = [C*, 2tt] and 
maps both I± onto the interval [E — W-i-,E — W-]. Under the hypotheses (Al) - (A3), the interval 
[E — E — W-] contains 2m + 2 branch points of the Bloch quasi-momentum, namely the points 
Ej for j = 2n — l,2n . . . 2(n + m). Therefore, on [0, 27r], one has 4m + 4 branch points of the complex 
momentum (p, j = 2n — 1, 2n, . . . 2(n + m), such that £ I± and £{Cf) = ^j- They satisfy the 
inequalities 

(5.2) < C2n-, <Gn<---< C^in+m) <C< Ct^n+m) < C^in+m)-l < ■ < Ctn-1 < 2^- 

The complex branch points all lie on Vl^~^(M). Reducing Y, at no loss of generality, we from now on 
assume that the strip Sy contains only the real branch points of the complex momentum. 

5.1.3. The sets Z and G. Consider the gaps and bands of the periodic operator (0.2). Let Z and G 
respectively be the pre-image (with respect to £) of the bands and the pre-image of the gaps. Clearly, 
ZuG = W~^{M.), and Z n G = 0. Connected components of Z and G are separated by branch points 
of the complex momentum. 

The set Z is 27r-periodic. Consider the part of Z situated in the inter- 
val [0, 27r]. Consider Z, the collection of subintervals of [0, 27r] defined in 
section 1.3.1. One has 



(5-3) = [C2j-i> C2j]> ij = [C2j> C2j-i]> j = 7i,n+l,...n + m. 

Under our assumptions on Y, the set Z f] Sy consists of the intervals of Z 
and their 27r-translates. 

Define the finite collection of intervals Q as in section 1.3.2. One has 
^^^^ 07 = (C2>C2"j+i): Qj = iCtj+vCtj)^ j = n,n + l,...n + m-l, 

0n-l = iCtn-l ~ 2'^'C2n-l)i Qn+m = (C2(n+m) ' ^2(n+m) ^ 

The set G Ci Sy consists of the following connected components: 

• the intervals g^, j = n, . . . n + m — 1; 

• the connected component containing ^ = (it is the union of the interval fln-i and the complex 
branch of VK~^(]R) n Sy passing through 0); 

• the connected component containing C = C* (it is the union of tli6 intcrvcil ^n-\-m 

and the 

complex branch of M^~^(]R) n 5y passing through ("*); 

• all the 27r-translates of the curves mentioned above. 

5.1.4. The main branch of the complex momentum. Introduce the main branch of the complex mo- 
mentum. In the strip {0 < Im(^ < Y}, consider the domain Dp between the complex branches of 
VF~^(R) beginning at and at C*. It is regular and £ = E — W{C) conformally maps it onto a domain 
in the upper half of the complex plane. We define the main branch of k by the formula 

(5.5) Kp{C) = kp{E-W{C)), C^Dp, 

where kp is the main branch of the Bloch quasi-momentum of the periodic operator (0.2) (see sec- 
tion 3.2). Clearly, ImKp > in Dp and the function Kp is continuous up to the boundary of Dp. Its 
behavior at the boundary of Dp reflects the behavior of kp along the real line. In particular, for each 
J = n, n -|- 1 . . . , n -|- m, it is monotonously increasing on and maps it onto [7r(j — 1), vrj]. 
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Figure 5: The branch 
points when m = 



5.1.5. The Stokes lines. Consider the Stokes lines beginning at C2n-i ^^'^ ^2n- '^p ^^^^ 
hn — [C2n-iiC2n]' ^^^^ interval is a Stokes line both for C,2n-i ^"^^ C2n- Pick one of these points. As 
W 7^ at this point, the angles between the Stokes lines beginning at it are equal to 2tt/3. So, one 
of the Stokes lines is going upwards, one is going downwards. These two Stokes lines are symmetric 
with respect to the real line. 

Denote by ai the Stokes line starting from C2^_^ going downwards, and denote by o"2 be the Stokes 
line beginning at and going upwards (see Fig. 6). 

The lines cri and a2 are vertical in Sy- Indeed, a Stokes line stays vertical as long as ImK ^ 0; the 
imaginary part of the complex momentum vanishes only on Z, and Z D Sy C M. 

Reducing Y if necessary, we can assume that cji and (T2 intersect the boundaries of Sy- 

5.2. Local construction of the solution /. We construct / on a local canonical domain. To 
construct a local canonical domain, we need a canonical line. To find a canonical line, we first build 
a pre-canonical line. 

5.2.1. Pre-canonical line. Consider the curve (3 which is the union of the Stokes lines cJi, [C2ri-i)C2ri] 
and (72. Let us construct a, a pre-canonical line close to the line (3. It goes around the branch points 
of the complex momentum as shown in Fig. 7. 

When speaking of Kp along a, we mean the branch of the complex momentum obtained of Kp by 
analytic continuation along this line (the analytic continuation can be done using formula (4.3)). 

Actually, the line a will be pre-canonical with respect to the branch of the complex momentum related 
to Kp by the formula: 

(5.6) Kp-n{n-l) if n is odd ^ 

I vrn — Hp it n is even 

In view of (4.5), k is indeed a branch of the complex momentum. We prove 

Lemma 5.1. Fix 5 > 0. In the 6 -neighborhood of (3, to the left of [3, there exists a, a line pre-canonical 
with respect to the branch k such that at its upper end IxaC, > Y, and at its lower end ImC < —Y. 

Proof. We consider only the case n odd. The analysis of the case of n even is similar. Note that, for 
n odd, formula (5.6) implies that K{C2n~i) = 0) and K(C2n) = So, the Stokes lines ai and a2 satisfy 
the equations Im ff- k = and Im r5_ (k — vr) = 0. 

f2n-l '>2n 

The pre-canonical line is constructed of lu and Id, two "elementary lines" that are segments of au and 
ad, two lines of Stokes type shown in Fig. 7. Let us describe them more precisely. 

Pick (^u, a point of the line Im^ = Y, to the left of /? close enough 
to it. The line Cu is the line of Stokes type Im {k — 7r)d( = Const 
containing (u- Note that all three curves, au, a2 and [C2n-i'C2nJ 
belong to the same family of curves Im j'^ [K—'K)dC, = Const. This implies 
that, if Cm is close enough to a2, then, below C,u, (^u goes arbitrarily close 
to £72 U [C2^_i,C2^] staying to the left of a2 and above [C2^_i,C2^]- We 
omit elementary details. 

Pick Cd, a point of the line Im Q = —Y, to the left of [3 close enough to it. 
The line ad is a line of Stokes type Im — 'K)dC, = Const containing 

Cd- Note that all the three curves ad, ai and al belong to the same family Figure 6: Stokes lines 
of curves Im ndC, = Const. This implies that, if Qd is close enough to 
(Ji, then, above C,d, o'd goes arbitrarily close to the line o"i U oT and stays 
to the left of it. We omit the details. 

Note that, by Lemma 4.2, the curves au and ad are tangent to the vector fields k — tt and k respectively. 
This implies in particular that each of the curves au and ad stays vertical as long as it does not intersect 
the set Z. 

If au and ad are chosen close enough to (3, they intersect one another in a neighborhood of C2n-i- 
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Indeed, consider first Uu- It goes from C,u downwards staying to the left of 02 and above [C,2ni^ C2n\- 
staying vertical, it has to intersect the Stokes line oT above C2n-i' beginning of oY- The intersection 
is transversal (as, first, oT is tangent to the vector field k, second, (T„ is tangent to the vector field 
K — TT, and, third, at the intersection point, ImK 7^ 0). If ad is sufficiently close to aTj then Gu also 
intersects ad transversally. 

We choose the intersection point as the lower end of lu and the upper end of Id- As au and ad are 
defined and are vertical somewhat outside S'y, we can assume that the upper end of lu is above the 
line Im^ = Y, that the lower end of Id is below the line ImC = —Y, and that both Id and lu are 
vertical. 

The lines lu and Id then form a pre-canonical line a; it can be chosen as close to the line (3 as desired 
and, in particular, inside the (5-neighborhood of /?. This completes the proof of Lemma 5.1. □ 



5.2.2. Local canonical domain and a solution f . By Proposition 4.1, as close to a as desired, one can 
find a canonical line 7. We construct 7 in the left part of the ^-neighborhood of (3. 

By Lemma 4.1, there \s K a canonical domain enclosing 7. We can assume that it is situated in 
arbitrarily small neighborhood of 7. So, we construct K in (the left part of) the (5-neighborhood of /?. 

By Theorem 4.1, on the canonical domain K, we construct /, a solution to (4.1) that has standard 
behavior / ~ exp J'' nd(^^ in K. We fix the normalization of this solution later. 

5.3. Asymptotics of / outside K. Recall that / is analytic in {Yi < Im^ < Y2}, the minimal strip 
containing K. 

5.3.1. The results. Denote by and the branch points situated on M and, respectively, closest to 
^2n-i °^ ^^'^ closest to °^ right. Let P be a regular domain obtained by cutting Sy 

along the Stokes lines ai, and a"2 and along the real intervals (— oo,Ca] and [Cb,oo), see Fig. 8. We 
prove 

Proposition 5.1. If 6 is sufficiently small, then, f has standard behavior 



in the whole domain T>. 

The rest of this section is devoted to the proof of Proposition 5.1. The proof is naturally divided into 
"elementary" steps. At each step, applying just one of the three continuation tools (i.e. the Rectangle 
Lemma, the Adjacent domain principle and the Stokes Lemma), we justify the standard behavior of 
/ on one more subdomain of D. Fig. 8 shows where we use each of the continuation principles. Full 
straight arrows indicate the use of the Rectangle Lemma, circular arrows indicate the use of the Stokes 
Lemma, and, in the hatched domains, we use the Adjacent Canonical Domain Principle. 

Again, the analysis of the cases n odd and n even are analogous. For sake of definiteness, we assume 
that n odd and consider only this case. We shall use 
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Figure 7: The construction of the pre-canonical curve 



(5.7) 



Lemma 5.2. If n is odd, then, in'D \ [C2n-i^C2n\' 

• to the left of the Stokes lines a\ and 02, one has Imn > 0; 

• to the right of these Stokes lines, one has Imn < 0. 

Proof. The sign of Iihk remains the same in any regular domain not intersecting Z. Moreover, the 
sign of Im k changes to the opposite one as C intersects a connected component of Z at a point where 
W / 0. So, in the connected subdomain of I? \ [C2n-i^C2n\ ^1 ^ [C2^-i'C2^] ^ '^2; one 

has luiK = Im(Kp — 7r(n — 1)) = Im/tp > (here, we have used (5.6) for n odd). To come from this 
subdomain to the connected subdomain of "D \ [C2^_i, C2^] to the right of cti U [C2^_i, C2^] U (T2, one has 
to intersect the interval [C2^_i, C2n] which is a connected component of Z. So, in the right subdomain, 
one has Im k < 0. □ 

5.3.2. Behavior of f between the lines 7 and (3. To justify the standard asymptotics of / in P between 
the lines 7 and (3, we use the Adjacent Canonical Domain Principle. Therefore, we need to describe 
a canonical domain enclosing 7 (more precisely, the part situated between 7 and /?). We do this by 
means of the Trapezium Lemma 4.4 (first statement). 

Let us describe the domain U and the curves 70, cTd and 
Gu needed to apply the Trapezium Lemma 4.4. 

The domain U . It is the domain bounded by /3, 7 and the lines 
ImC, = Const containing the ends of 7. 

In view of Lemma 5.2, choosing the ends of 7 closer to the lines 
ImC^ = ity if necessary, we can assume that ImK > in the 
domain U . 

The line Cu- As the line Cu, we take the line which belongs to 
the family Im ndQ = Const and intersects (3 at C,u satisfying 
Im C,u = Y . Recall that 7 is constructed in the ^-neighborhood 
of (3 where 5 can be fixed arbitrarily small. One has 

Lemma 5.3. The line cj„ enters U at Cu o,nd goes upwards. If 5 is sufficiently small, then, intersects 
7 at an internal point of 7. 

Proof. The main tool of this proof is Lemma 4.2. Below, we use it without further notice. Recall that, 
above the real line, (3 coincides with the Stokes line a2. So, it is tangent to the vector field k(^) — vr. 
The line ad is tangent to the vector field k(C)- In U, and in particular at Qu, one has ImK > 0. 
Therefore, at C,u, the tangent vector to (3 (oriented upwards) is directed to the right with respect to 
the tangent vector to Uu (oriented upwards). So, Uu enters U at Cm going upwards. As Im/t / in [/, 
it stays vertical in U . Note that au is independent of 6. So, if 5 is small enough, then au intersects 
7. □ 

The line cr^. It is the line of Stokes type Im j'^ KdQ = Const that intersects [3 at Q satisfying 
Im^d = —Y- One has 

Lemma 5.4. The line cj^ enters U at Cd cmd goes downwards. If 5 is sufficiently small, then 
intersects 7 at an internal point. 

The proof of this lemma being similar to the proof of Lemma 5.3, we omit it. 

The line 70. We choose 6 so that both au and ad intersect 7. Then, 70 is the segment of 7 between 
its intersection points with ad and cr^. 

5.3.3. Describing the curve 7. Let us describe the canonical line 7 needed to apply the first variant of 
the Trapezium Lemma. As for 7, using Lemma 5.1, we can construct 7 so that it be arbitrarily close 
to P and strictly between 70 and /?. 

As au and ad intersect 7 and /3, they also intersect 7. 
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Figure 8: How to "continue" the 
asymptotics of / 



5.3.4. Completing the analysis. By the Trapezium Lemma, the domain bounded by 70, cr^, Cd and 7 
is a part of the canonical domain enclosing 79. So, by the Adjacent Canonical Domain Principle, / 
has the standard behavior here. 

As 7 can be constructed arbitrarily close to (3, we conclude that / has the standard behavior in the 
whole domain bounded by /3, 7 and the lines |Im(^| = Y. 

5.4. Behavior of / to the left of 7. We justify the standard behavior of / in P to the left of 7 by 
means of the Rectangle Lemma. 

Let us describe R, the rectangle used to apply Lemma 4.3: it is the part of V between the canonical 
line 7i = 7 and a vertical line in D, say 72, staying to the left of 7 and going from Im^ = —Y to 
lmC = Y. 

By Lemma 5.2, in R, the imaginary part of k is positive. Moreover, / has standard behavior in a 
neighborhood of 7. So, by the Rectangle Lemma, / has standard behavior in R. 

Pick £ T> to the left of 7. As the vertical curve 72 C P can be taken so that C, G R, we see that / 
has standard behavior in the whole part of T> situated to the left of 7. 

5.5. Behavior of / to the right of ui. First, by means of the Stokes Lemma, Lemma 4.5, we show 
that / has standard behavior to the right of ai in its small neighborhood. 

Let Vi be a sufficiently small constant neighborhood of cji and set Vi = Vi D Sy- Show that / has 
standard behavior in Vi\ai. The Stokes lines di, oT and [C2n-ii C2n\ divide Vi into three sectors. By 
the previous steps, we know that / has the standard behavior in Vi outside the sector S bounded by 
the Stokes lines ai and [C2n-i'C2n]- The Stokes line ai is vertical. By Lemma 5.2, in Vi, to the left 
of o"!, the imaginary part of the complex momentum is positive; thus, the Stokes Lemma implies that 
/ has standard behavior inside Vi \ cJi . Recall that the leading term of the asymptotics of / in 5 is 
obtained by analytic continuation from the rest of Vi around the branch point C2n-i avoiding ai. 

Having justified the standard behavior of / to the right cJi in a small neighborhood of di, we justify 
it in in the rest of the subdomain of D n { — Y < lm(^ < 0} situated to the right of ai by means of 
the Rectangle Lemma. The argument is similar to that carried out in subsection 5.4. So, we omit the 
details. 

5.6. Behavior of / along the interval (C2n-i'C2n)* previous steps, we have justified the 
standard behavior of / both above and below the real interval (C2ri-i' C2ri)- show now that / has 
standard behavior also in this interval. 

By section 5.5, we know that / has standard behavior in a neighborhood of C2n-i along o"i. 
Moreover, / can not have the standard behavior in a neighborhood (as ^2^ is a branch point). 
Hence, there exists a E (C2ri-i'C2ri] such that / has the standard behavior in a neighborhood of any 
point in (C2n-i''^)- Assume that a < C,2n- a be a segment of the line ReC = a connecting a point 
ai S C_ to a point 02 € C+. One has < K(a) < vr. This implies that, if a is small enough, then, 
a is a canonical line. The solution / has the standard behavior to the left of a. By the Adjacent 
Canonical Domain Principle, / has standard behavior in any local canonical domain enclosing a, thus, 
in a constant neighborhood of a. So, we obtain a contradiction, and, thus, a = (^^. This completes 
the analysis of / along (C2n~i'C2n)- 

5.7. Behavior of / to the right of c72. One studies / to the right of (T2 in the same way as we have 
studied it to the right of ai: first, using the Stokes Lemma, one justifies the standard behavior to the 
right of £72, in a small neighborhood of CJ2, and, then, applying the Rectangle Lemma, one proves that 
/ has the standard behavior in the rest of the subdomain of T) to the right of this neighborhood. We 
omit further details. 

The analysis of / to the right of (72 completes the proof of the Proposition 5.1. □ 

5.8. Normalization of /. To fix the normalization of the leading term of the asymptotics of /, we 
choose the normalization point Co in ^ and, in a neighborhood of we choose a branch of the function 

19 



^/k'{£{C)) in the definition of 

As the normalization point, we take Co such that 



(5-8) C2n-1 < Co < Can 

Inside any spectral band of the periodic operator, k'{£) does not vanish, and there are no poles of the 
Bloch solution iIj{x,£). So, £{C,o) ^ P U Q, and the solution "if^ is well defined. 

To fix the branch of ^/W, we note that, inside any spectral band of the periodic operator, the main 
branch of the Bloch quasi-momentum, kp, is real and satisfies k'^ > 0. So, we fix -v/a7 so that 



(5.9) yW{0)>0, C2n-l<C<C2n- 

6. The consistent basis 

Up to now, we have constructed /, one consistent solution to (4.1) with known asymptotic 
behavior in the domain D. We now construct another consistent solution /* so that {f,f*) form a 
consistent basis. 

6.1. Preliminaries. For C G T^*, the symmetric to D with respect to the real line, we define 



(6.1) nx,C,E)=f{x,C,E). 

As W is real analytic, the function /* is also a solution of (4.1); it satisfies the consistency condition 
as / does. In the next subsections, we first study its asymptotic behavior; then, we compute the 
Wronskian w{f, /*). We show that, in Sy, it has the form Const (1 + o(l)). Here, Const is a non-zero 
constant, and o(l) is a function which can depend on C- Finally, we modify the solution / so that it 
still have the standard behavior in V, and w{f, /*) be constant. 

6.2. The asymptotics of /*. Note that VnV* contains the interval 3 = {C2n-i^ C^n) ^- One has 
Lemma 6.1. In D* , the solution f* has the standard behavior 

(6.2) /~e-^^Co''*'^^^_,,(x,C,Co), 
where 

• is the branch of the complex momentum which coincides with k. on ^ and is analytic in T>* , 

• is the canonical Bloch solution which coincides with (corresponding to from the 
asymptotics of f) on i and is analytic in D* . 

Proof. Recall that, by Proposition 5.1, / has the standard behavior (5.7) in the domain V. 
The statement of Lemma 6.1 follows from Proposition 5.1, the definition of /* and the relation 



(6.3) exp KdCj ^+ix, C, Co) = exp (^-'- K,dC^ C, Co). C G 

Let us prove this relation. As both the right and left hand sides of (6.3) are analytic in C, it suffices 
to check (6.3) along the interval 3. Recall that the interval [Can-i'Can] ^ connected component of 
Z. This implies that 



(6.4) ^(0 = ^^(0, i^+{x,£{0) = i^-{x, £{()), Cea, 

where ip±{x,£) are two different branches of the Bloch solution ip{x,£). 

As Co satisfies (5.8), relation (6.3) follows from the first relation in (6.4) and the relation 



(6.5) ^'+(x,C,Co) = ^-(a^,C,Co), CG3. 

To check (6.5), we recall that ^± are defined in section 4.1.2 by formula (4.7). Therefore, relation (6.5) 
follows from (5.9), the second relation in (6.4) and the last property of uj listed in the section 3.3. This 
completes the proof of Lemma 6.1. □ 
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6.3. The Wronskian of / and /*. The solutions / and /* are analytic in the strip Sy- Here, we 
study their Wronskian. As both / and /* satisfy condition (4.2), the Wronskian is e-periodic in 
One has 

Lemma 6.2. The Wronskian of f and f* admits the asymptotic representation: 

(6.6) w{f,n=w{^+,-i>-)\c=Co+9, C^Sy. 

Here, g is a function analytic in Sy, such that, along the real line, Reg' = 0. Moreover, g = o(l) 
locally uniformly in any compact of Sy provided that E is in a sufficiently small complex neighborhood 
ofEo. 

Remark 6.1. Note that 

(1) u;(^'+,^'_)|f=<;o / (as £:(Co) -PU Q, see (5.8) and the comments to it); 

(2) ?i;(^'+,^'_)|(^=^o e iM (due to (6.5)). 

Proof. The domain D f] D* contains the "rectangle" R bounded by the lines ai U oi, (T2 U 02 and 
Im(^ = zizY. So, in R, the solutions / and /* have the standard behavior (5.7) and (6.2). Consider 
the functions k^, and from (6.2). Their definitions, see Lemma 6.1, imply that 

= K, ^'_,* = ^'_, C £ R. 

Using this information and (5.7) and (6.2), one obtains 

(6.7) w{f,f*)=w{^+{;C),^^{;0)+9, 9 = o{l), ( ^ R- 

Being obtained using standard behavior, this estimate is uniform in in any compact of R provided 
£J be in a sufficiently small neighborhood of Eq. By (4.8), the first term in the left hand side of (6.7) 
coincides with the first term in (6.6). So, we only have to check that g has all the properties announced 
in Lemma 6.2. As w{f,f*) is e-periodic, so is g. Furthermore, ig is real analytic as iw{f,f*) and 
iw{^-^, ^'_)|^=^Q are. This completes the proof of Lemma 6.2. □ 

6.4. Modifying /. As g, the error term in (6.6) may depend on C, we redefine the solution /: 

/ := f/u, V = yr+^M^T^^iy^. 

In terms of this new solution /, we define the new /* by (6.1). These are the basis solutions the 
monodromy matrix of which we shall study. For these "new" functions / and /*, we have 

Theorem 6.1. The solutions f and f* satisfy the condition (4.2) and 

(6.8) «;(/,/*) =^*+,*_)|c=^,. 

Moreover, f has the standard behavior (5.7) in T>, and f* has the standard behavior (6.2) in T>* . 

Proof. Let C be in a fixed strip {yi < ImC < 2/2} C Sy, and let e be sufficiently small. We use 
Lemma 6.2 and Remark 6.1. Recall that g is e-periodic in So, 1^ is e-periodic, and / and /* remain 
consistent. Furthermore, note that is real analytic. This implies (6.8). Finally, as = 1 -|- o(l), the 
new solutions / and /* still have the "old" standard asymptotic behavior in V and V* respectively. □ 



7. General properties of the monodromy matrix for the basis {/,/*} 

In the previous section (see Theorem 6.1), we have constructed (/, /*), a consistent basis of solutions 
of (4.1). If we return to the variables of the initial equation (0.1), we get a consistent basis of (0.1). 
The matrix discussed in Theorem 1.4 is the monodromy matrix obtained for this basis. In this short 
section, we check some of its properties. 

Instead of coming back to the initial variables, we continue to work in the variables (5.1). The 
definition of the monodromy matrix (1.6) takes the form 

(7.1) F(x,C + 27r)=M(C)F(x,C), F = (^jj^^^^^y 
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and the matrix M becomes e-periodic in ^. 

As the basis solutions / and /* are related by (6.1), the monodromy matrix has the form (1.10). 
The definition of the monodromy matrix (7.1) implies that 

(7 9\ fn-M (r^ ^(/(^ + 2vr,C),r(x,C)) yn-M (n ^(/(^,C),/(^,C + 2^)) 

Finally, we note that the monodromy matrix is analytic in in the strip Sy and in in a constant 
neighborhood of £"0. Indeed, as the solutions / and /* are analytic functions of both variables, so 
are the Wronskians in (7.2). Moreover, by (6.8), the Wronskian in the denominators in (7.2) does not 
vanish. Hence, we have proved the 

Lemma 7.1. The monodromy matrix corresponding to the basis {/,/*} satisfies (7.1), is analytic and 
e-periodic in ( £ Sy, analytic in E in a constant neighborhood of Eq and has the form (1.10). Its 
coefficients are given by (7.2). 



8. General asymptotic formulas 

To compute the asymptotics of the monodromy matrix defined above, we only need to compute the 
Wronskians in the numerators in (7.2). These Wronskians depend on C and have different asymptotics 
in the lower and upper half planes. Rather than repeating similar computations many times, in the 
present section, we obtain a general asymptotic formula for the Wronskian of two solutions having 
standard behavior. 

8.1. General setting. In this subsection, we do not suppose that W be periodic. Fix E = Eq. 
Assume that h and g are two solutions of (4.1) having the standard asymptotic behavior in regular 
domains Dh and Dg-. 

(8.1) /i~e^4"''''^M/^(x,C), 5~e^^c«"«^^vf9(x,C). 

Here, Kh and Kg are branches of the complex momentum analytic in Dh and Dg, and ^g are 
canonical Bloch solutions defined on Dh and Dg, and Cfe and are the normalization points for 
h and g. 

As the solutions h and g satisfy the consistency condition, their Wronskian is e-periodic in We now 
describe the asymptotics of this Wronskian and of its Fourier coefficients. We first introduce several 
simple useful objects. 

Below, we assume that Dg n Dh contains a simply connected domain, say d. 

8.1.1. Arcs. Let 7 be a regular curve going from (^g to C/, in the following way: staying in Dg, it goes 
from (g to some point in d, then, staying in Dh, it goes to Ch,. We say that 7 is an arc associated to 
the triple h, g and d. 

As d is simply connected, all the arcs associated to one and the same triple can naturally be considered 
as equivalent; we denote them by ^{h,g,d). 

Let us continue Hh and Kg analytically along ^[h,g,d). The analysis performed in section 4.1.1, 
see (4.5), yields, that, for V a small neighborhood of 7, one has 

(8.2) Kg{C) = (TKh{C,) + 2TTm, m£Z, cjG{-1,+1}, C^V. 
We call a = a{h, g, d) the signature of 7, and m = m{h, g, d) the index of 7(/i, g, d). 



8.1.2. Meeting domains. Let d be as above. We call d a meeting domain, if, in d, the functions ImK/j 
and Im Kg do not vanish and are of opposite sign. 

Note that, for small values of e, the increasing and decreasing of h and g is determined by the 
exponential factors ^ '^'^^ . So, roughly, in a meeting domain, along the lines ImC = Const, the 
solutions h and g increase in opposite directions. 
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8.1.3. The amplitude and the action of an arc. We call the integral 

S{h,g,d) = / KgdC 

the action of the arc 7(/i, g, d). Clearly, the action takes the same value for equivalent arcs. 

Assume that £{(,) ^ P^Q along 7 = 7(^1, g^ d). Consider the function qg = \/k'{£{C)) and the 1-form 
Q.g{£{C,)) in the definition of ^g. Continue them analytically along 7. Put 

(8.3) A(/i,<7,7) = {%l<lh\=CH 

We call A the amplitude of the arc 7. The first three properties of Vt listed in section 3.3 imply 
Lemma 8.1. The amplitudes of two equivalent arcs "f{h,g.,d) coincide. 

8.1.4. Fourier coefficients. Let S[d) be the smallest strip of the form {Ci < IvaC, < C2} containing 
the domain d. One has 

Proposition 8.1. Let d = d{h,g) be a meeting domain for h and g, and m = m{h,g,d) be the 
corresponding index. Then, 

(8.4) w{h,g) = Wme"-^^<~'^^\l + o{l)), C G S{d), 
and Wm is the constant given by 

(8.5) Wm = A{h, g, d) el ^('^'^^''^) w{-^+{; (h), *-(•, Ch)), 

where = Oind ^_ is "complementary" to The asymptotics (8.4) is uniform in Q and E 

when Q stays in a fixed compact of S{d) and E in a small enough constant neighborhood of Eq. 

Note that the factor Wm is the leading term of the asymptotics of the m-th Fourier coefficient of 
w{h,g). 

Proof. For E Dg, let 7g(C) be a curve in Dg from (g to Similarly, define 7h(C)- 
First, we check that, for C ^ d,, one has 

(8.6) '^^ = ^'-^'^''^^ e^ (e~^ ^'"^('^^ '^'^^ , 

(8.7) ^g{x,C)=A{h,g,d)-i'^ix,C), 

where is the canonical Bloch solution "complementary" to ^+ = in a neighborhood of 7/i(C)- 
As d is a meeting domain, in a neighborhood of 7(/i, g, d), one has 

(8.8) Kg = —Kh + 2'Km 
This implies relation (8.6). 

Check (8.7). Let 7 = ^{h,g,d) be an arc such that S{'j) H {P U Q) = 0. Continue qg, Qg and ^pg 
analytically along the arc 'j{h,g,d). Note that qg and q^ are two different branches of the function 
^k'(E — W((y). So, they differ at most by a constant factor. Therefore, in a neighborhood of Qh-, we 
get 

(8.9) ^g{xX) = A{h,gA)qh{e{Q)e^'^''''''''^g{x,£{Q). 

Now, recall that, in a neighborhood of C^, there are only two branches of Q and ■(/''• Denote them 
by and Vl± so that -0+ = V'/i a-^^d = fi/^. Then, either '^g = ip- and Q.g = or ipg = ip^ 
and Qg = To choose between these two variants, we recall that the Bloch quasi-momentum of 
a Bloch solution is defined modulo 27r. Note that Kh is the Bloch quasi-momentum of ip^, and Kg 
is the Bloch quasi-momentum of ipg. By (8.8), we get Kg = — k^, mod27r. So, Kg must be the Bloch 
quasi-momentum of ip-. Thus, we have ipg = tp- and Qg = and (8.9) implies relation (8.7) in a 
neighborhood of Q^. By analyticity, it stays valid in d. 

As d C Dh n Dg, both h and g have standard behavior in d. Substituting the asymptotics of / and g 
into w{f,g), and using (8.6) and (8.7), one easily obtains 

(8.10) w{h, g) = A{h, g, d) el ^('^'f''^) [w{^+{; C), ^-(•, 0) + «(!)), C e d. 
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As w{'^+{-, C), ^-(-j C)) is independent of ( and e and is non-zero (see (4.8) and comments to it), we 
get (8.4). As this asymptotic was obtained using the standard behavior of h and g, it has all the 
announced uniformity properties. This completes the proof of Proposition 8.1. □ 

8.2. The index m and the periods when W is periodic. Here, we only assume that W is 27r- 
periodic and real analytic in (i.e. we do not assume anything on the critical points of W), and that 
E is fixed. We describe the computation of the index m in the special case that one encounters when 
computing monodromy matrices. 

8.2.1. Periods. Pick a regular point. Consider a regular curve 7 going from (q to Co + 2vr. Fix k, 
a branch of the complex momentum that is continuous on 7. We call the couple (7, k) a period. 
Let (71, Ki) and (72, K2) be two periods. Assume that one can continuously deform 71 into 72 without 
intersecting any branching point. This defines an analytic continuation of ki to 72. If the analytic 
continuation coincides with K2, we say that the periods are equivalent. 

Consider the branch k along the curve 7 of a period (7, k). In a neighborhood of Co) the starting point 
7, one has 

(8.11) k(C + 27r) = fJK(C) + 27rm, cj G {±1}, m £ Z. 

The numbers a = cr{'y, k) and m = m{^, k) are called the signature and the index of the period (7, k). 
The numbers m (resp. a) coincide for equivalent periods. 

Recall that G is the pre-image with respect to £ of the spectral gaps of Hq. One has 

Lemma 8.2. Let (7, k) he a period such that 7 begins at a point Co G- Assume that 7 intersects 
G exactly N times (N G W ) and that, at all intersection points, W / 0. Let ri, r2, . . . , r^ he the 
values that Rck takes consecutively at these intersection points as C, moves along 7 from C,q to C,q + 2tt. 
Then, 

(8.12) ,7(7, k) = (-1)^, m(7,K) = -(r^-r^_i+r^_2---- + (-l)^-Vi). 

vr 

Proof. The image f (7) of 7 by <S : ( E — W{C) is a closed curve that starts and ends at Sq = £{Co)- 
We consider the curve as open at £q. Along 7, we can write k{() = k{E — W{C,)) where /c is a 
fixed analytic branch of the quasi-momentum. So, h{Cq) and k(Co + ^vr), the values of the complex 
momentum at the ends of 70, are related by the same formula as k^, and fcg, the values of k at the 
beginning and the end of curve ^(70)- 

Since W' / at the points of intersection of 70 and G, intersects exactly N times spectral gaps 

of the periodic operator. 

As the values for both m and a coincide for equivalent periods, it suffices to construct Qq so that 
Im£:o / 0. 

Assume that a continuous curve begins at £.q, goes along a strait line to one of the ends of a gap, then 
goes around this gap end along an infinitesimally small circle, and returns back to £q along the same 
strait line. We call such a curve a simple loop. We distinguish the end and the beginning of the loop 
considering it as open at its endpoints. As ImS-Q 7^ 0, any simple loop intersects only one gap, namely, 
the gap around the end of which it goes. 

Recall that the ends of the gaps coincide with the branching points of the Bloch quasi-momentum, 
and that these branching points are of square root type. So, in a neighborhood of a branching point, 
the corresponding branches of the Bloch quasi-momentum satisfy the relation 

(8.13) ki{E) + k2{E) = 2r, 

where r is the common value of these branches at the branching point. Note that r is equal to the 
value of the real part of any of these branches on the spectral gap beginning at the branching point. 
On a simple loop, fix a continuous branch of the quasi-momentum. Clearly, formula (8.13) also relates 
the values of the quasi-momentum at the ends of the loop when r is the value of the quasi-momentum 
at the branching point inside the loop. 

Recall that k can be analytically continued onto the whole complex plane cut along the spectral 
gaps of Hq. Therefore, the value of k at the end of <f (70) is equal to the value of k at the end of 
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the curve consisting of N simple loops and going successively around the branch points of k with 
k = ri,r2, ■ ■ - vn. So, taking (8.13) into account, we get 

ke = (-l)^^b + 2(rAr - rAr„i + r^Tj - . . . Ti) . 
This implies (8.12) and completes the proof of Lemma 8.2. □ 

8.2.2. Indices of periods. Let us come back to the computation of the index m{h,g,d). One has 
Lemma 8.3. Let 7 = 7(/i, g, d) he an arc such that C,h = Cg + 2vr. //, in a neighborhood of Qg, 

(8.14) Kg{C) = s-Kh{C + 2n), 
where s is either "+" or then, 

(8.15) a{h,g,d) = s ■ a{'j,Kg), m{h,g,d) = m{'j,Kg). 

Proof. The pair {j{h, g,d), Kg) is a period. So, in a neighborhood of Cg, one has Hg{C + 2vr) = 
a{'y, K,g)K,g{() + 2iTm{'-),Hg). This and (8.14) imply that Hg{C) = s a{'y, Kg)Kh{C) + 2iTm{'j, Kg) in a 
neighborhood of (h- This implies the relations (8.15). □ 

9. ASYMPTOTICS OF THE MONODROMY MATRIX 

We now compute the asymptotics of the coefficients a and b of the monodromy matrix for the 
basis {/,/*}; in particular, we prove formulae (1.11) and (1.12). We concentrate on the case n odd. 
The computations for n even being similar, we omit them. 

Recall that a and b are expressed via the Wronskians by formulae (7.2). We compute these Wronskians 
(and, thus, a and b) using the construction from section 8. 

9.1. The asymptotics of the coefficient b. By (7.2), we have to compute u)(/(-, C), /(•, C + 2vr)). 
One applies the constructions of section 8. Now, one has 



(9.1) h{x,C) = f{x,C), gix,C) = {Tf){x,C) where (r/)(x, C) = /(x, C + 27r); 

(9.2) Dh = V, Dg = V- 2tt; 

(9.3) Ch = Co, = Co - 27r; 

(9.4) nh{0=K{C), Kg{C)=K{C + 2^). 



9.1.1. The asymptotics in the strip < ImC < Y- Let us describe do, the meeting domain, and 
7o = l{f f ,dQ), the arcs used to compute w{f,Tf) in the strip {0 < ImC < 

The meeting domain, do is the subdomain of the strip < Im C < ^ between the Stokes lines (T2 — 27r 
and fT2. Indeed, it follows from Lemma 5.2 and (9.4) that, in do, one has ImK^ = —luiKh < 0. 
The arc. 70 connects the point C3 to C/i- In view of (9.3), it defines the period (70, Kg). 
The index m{f,Tf,do). In view of (9.4), the arc 70 satisfies the assumption of Lemma 8.3. So, 
m{f,Tf,dQ) = m{'yo,Kg). Due to (9.4), m{f,Tf,do) = m(7o + 2'k,k). To compute this integer, we 
use Lemma 8.2. Therefore, we have to compute k at the intersections of 70 + 27r and G, the pre-image 
of the spectral gaps of Hq. The set G Ci Sy is described in section 5.1.3 where we have listed all its 
connected components. 

Recall that m takes the same value for all the periods equivalent to (70 + 2tt,k). We can deform 
7o + 27r into a curve, say ■j C V, so that 

• 7 be to the left of the complex branch of VF~^(M) starting at 2tt and staying in the upper 
half-plane, 

• (7, k) be a period equivalent to (70 + 2-ir, k), 

• 7 have the following intersections with the connected components of G (for m = 0, this curve 
is shown in Fig. 9): it once intersects the complex branch of VK~^(M) going upwards from 
0, once the interval (C2^,Cfe) (t^ie point Cfe is defined in section 5.3.1), and, once the complex 
branch of VF~^(]R) going upwards from C*- 
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Figure 9: Curves equivalent to 70 + 2tt and 71 + 27r (when m = 0) 



Recall that ReK(^) is constant on any connected component of G. Therefore, 

m(7o + 27r,K) = m{j,K) = ^ (ReK(C)|^-^_^ - ReK(C)|^^- + Re K(C)|^g0„^^+io) 

= i(o-vr+ReK(C)|^,,„^^+,o)- 

To compute the last term in this formula, we recall that, in -Dp, left, the part of the domain Dp (see 
section 5.1.4) situated to the left of the Stokes line (T2, one has k = Kp — 7r(n — 1). In the domain 
-Dp,right ) the part of Dp situated to the right of (J2 , k is obtained by analytic continuation from Dp jeft 
around the branch point passing below this branch point. As Kp((^^) = vrn, for ( £ -Dp^ight, one 
has k{() = {2-Kn — Kp{()) — 7r(n — 1) = 7r(n + 1) — Kp{(). Along Qn+m, one has RcKp = 7r(n + m); 
hence, we get 

(9.5) m(7o + 27r, k) = -(0 - vr + [7r(n + 1) - 7r(n + m)]) = -m. 

vr 

The result. Now, Proposition 8.1, formula (7.2) for b and formula (6.8) imply formula (1-11) for b with 

(9.6) = A{f, Tf, do) el^(/''^/''^o)+^, {Tf){x, () = f{x, C + 2n). 

9.1.2. Asymptotics of b below the real line. Describe di, the meeting domain, and compute the index 
m{f,Tf,di) needed to get the asymptotics of w{f,Tf) in the strip {—Y <lm( < 0}. 
The meeting domain, di is the subdomain of the strip —Y < Im C < situated between the Stokes 
lines o"! — 27r and ai. 

The index. The arc 71 = 7(/, T/, di) again defines a period, and m{f,Tf,di) = m{'ji + 2tt,k). The 
curve defining a period equivalent to (7i+27r, k) is shown in Fig. 9. As in the sequel of this computation 
we only use this curve, we call it 71 + 2tt. To compute the index of this period, we compute Re k at 
the intersections 71 + 27r and G. 
We can assume that 71 + 27r satisfies: 

• it is situated to the left of the complex branch of VF~^(R) in C_ starting at 27r, 

• it has the following intersections with G: it once intersects the complex branch of Vl^~^(M) 
going downward from 0, once the interval fln-i and once the complex branch of 1^~^(M) going 
downward from C,* . 

We get 

m(7i +27r,K) = ^ . 

We have used the fact that the interval Qn-i and the complex branch of VF~"'^(IR) going downward 
from belong to the same connected component of VF~^(]R). To finish the computation, we introduce 
the domain the symmetric of Dp with respect to the real line. In -D* ^.jgj^^, the part of this domain 
situated to the right of cji, k can be viewed as the analytic continuation of Hp — 7r(n — 1) from Dp 
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across the interval 3„ . Along the interval 3„ , Kp is real. So, for C, G -D* j-jg^tJ '^(C) = ^p(C) ~ '''"('^ ~ 1)) 
and 

(9.7) m(7i +27r,K) = ^ (^Re Kp(C)|(.gg^^^_^iO - vr(n - 1)) = :^(7r(n + m) - 7r(n - 1)) = m + 1 
The result. Now, Proposition 8.1, formula (7.2) for h and (6.8) imply formula (1.12) for b with 

(9.8) Wi =^(/,T/,di)e^^(^™) (r/)(x,C) =/(x,C + 2vr). 
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Figure 10: Curves 70 + 27r and 71 + 27r (when m = 0) 



9.2. The asymptotics of the coefficient a. The computations of the coefficient a following the 

same scheme as those of 6, we only outline them. Now, 

(9.9) h = r, g = Tf; Dh=V*, Dg =V-27r; 

(9.10) Ch = Co, Cg = Co - 27r; 

(9.11) KhiC) = -K{0, C^Dh, AJg(C) = At(C + 2vr), QCiDg. 

Recall that the complex momentum is real on 3", and that Co £ In- This and relations (9.11) imply 
that 

(9.12) ^^(^) = _^^(^ + 2^), C~Cr 

9.2.1. The asymptotics of a above the real line. In this case, do, the meeting domain, is the subdomain 
of the strip {0 < Im C < situated between the lines o"2 — 2tt and oT (which is symmetric to o"i with 
respect to M). The arc ^{f* ,Tf,do) defines a period (70, Kg); the curve 70 + 27r is shown in Fig. 10. 
In view of (9.12), one is again in the case of Lemma 8.14, and, by means of Lemma 8.2, one obtains 
m(/*,r/, do) = ^ilo + ^vr, k) = —m. This yields formula (1.11) for a with 

(9.13) a_„ = Air,Tf, do) Sir,Tf,<io)+^-^ _ 

9.2.2. The asymptotics of a below the real line. In this case, di, the meeting domain, is the subdomain 
of the strip {—Y < ImC < 0} situated between the lines 02 (symmetric to cj2 with respect to M) and 
0"! — 27r. The arc ^{h,g,di) defines a period (71, k^); the curve 71 + 27r is shown in Fig. 10. One 
obtains m{f*,Tf,di) = m + 1. This yields formula (1-12) for a with 



(9.14) a^+i = A{f*,Tf,di)e 
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10. ISO-ENERGY CURVE 



The iso-energy curve F is defined by (0.4). In this formula, £{■) is the dispersion law for the periodic 
operator (0.2) i.e. the function inverse to the Bloch quasi-momentum (£ = £{k) if and only if k is the 
value of one of the branches of k when the spectral parameter is equal to £). We begin with a simple 
general observation: 

Lemma 10.1. The iso-energy curve T is In-periodic in C- and k- directions; it is symmetric with 
respect to any of the lines k = irm, ?n G Z. 

Proof. The periodicity in follows from the one of W. Fix ko G C. The list (3.2) shows that £{k) 
takes the same value for all k = ak^ + 27rm, where a G {il}, and m G Z. This implies the periodicity 
and the symmetries in k. □ 

Now, for W satisfying (H) and for E in J, an interval satisfying (Al) - (A3), we discuss the iso-energy 
curves (0.3) and (0.4) and obtain the estimates (1.13). 

10.1. Real iso-energy curve: the proof of Lemma 1.1. A point {(, k) G belongs to Fk if and 
only if K is the value of one of the branches of the complex momentum at Recall that the intervals 
I ^ Z are pre-images (with respect to of spectral bands. The complement of these intervals in 
(0, 27r) is mapped by E into spectral gaps. So, on (0, 27r), k((^) takes real values only on the intervals 
of Z. Therefore, in the strip {0 < k < 27r}, the connected components of Fk are situated above the 
intervals I & Z ("above" refers to the projection 11 : (C, k) G ^ C G M). 

Pick J G {n, n + 1 . . . n + m}, and a G {±}. Consider the part of F]r above the interval i := 3J. 
Recall that £ bijectively maps 3 onto the j-th spectral band of Hq. So, there exists kq, a branch of 
the complex momentum, continuous on the interval 3, and mapping it monotonously onto the interval 
[7r(j - l),vrj] so that ft:o(C2i-i) = ^(j - 1) and Ko(C2i) = 

On the interval [7r(j — l),7rj], let be the inverse of ^ 1— > k;o(C)- We continue k to the real 

line so that it be 2tt periodic and even. Recall that all the values of all the branches of the complex 
momentum at C S 3 are given by the list (4.5). This and the definition of Z^ imply that the points of 
F]R above 3 are points of the graph of and reciprocally. 

All the properties of the function announced in Lemma 1.1 follow directly from this construction. 

To prove that the connected components of Fk depend continuously on E, it suffices to check that 
each of the functions Z^ depends continuously on E £ J. Pick ^ £ Z. As W'{C) 7^ for all C ^ h the 
continuity oi E Z^ immediately follows from the Local Inversion Theorem and the definition of the 
iso-energy curve (0.3). This completes the proof of Lemma 1.1. □ 

10.2. Loops on the complex iso-energy curve. Here, we discuss closed curves in F. 

10.2.1. An observation. We define the intervals g G ^ as in section 1.3.2. We shall use 

Lemma 10.2. Pick Q G G. Let V{q) be complex neighborhood of q sufficiently small so that it contains 
only two branch points of k, namely, the ends ofQ. Let k. be a branch of the complex momentum analytic 
in a sufficiently small neighborhood of a point of V{q) \ q. Then, k can be analytically continued to 
the domain V{q)\q to a single valued function. The analytic continuation satisfies 

(10.1) ^(0=^, C^V{q)\q. 

Proof. We can continue k to a branch of the complex momentum analytic in V'{q) the simply 
connected domain obtained from V{q) by cutting it, say, along M from the right end of g to -|-oo. It 
suffices to check that the values of k at the edges of the cut coincide. 

The set (M H V{q)) \ q consists of two intervals. Each of them belongs to Z (the pre- image of the 
spectral bands with respect to £). So, k is real both on the left of these two intervals and at the edges 
of the cut. As K is real on the left interval, one has (10.1) in V' . So, the values of k, on the edges of the 
cut satisfy k{C, + «0) = k(^ — iO), and, therefore, being real, coincide. This implies Lemma 10.2. □ 
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10.2.2. The loops. Pick Q £ Q. On V{q) \ g, fix acq, a single valued analytic branch of the complex 
momentum. Consider G{q) C V{g) \ 0, a curve going once around the interval g. One has 

Lemma 10.3. For each a G {±1} and m G Z, the curve 

G(m,a){9) = {{(,1^) ■■ K = o-Ko(C) + 27rm, C G G(g)}, m G Z, 

is a closed curve onV. It connects the two connected components ofT^ that project onto the intervals 
of Z DM. adjacent to q . 

Proof. As kq is univalent on G(g), the curve G{ofl) is a closed curve on T. This and Lemma 10.1 
imply that all the curves are loops in F. As G{q) intersects the intervals of Z n M adjacent to 

g, G connects the two connected components of Fr that project onto these intervals. □ 

10.3. Tunneling coefficients. Pick Q £ Q. Fix an analytic branch k, of the complex momentum on 
V{q). Define the action S{q) = i /(^(g) ^f^C- To study its properties, we use 

Lemma 10.4. Let E £ J. If G{q) is positively oriented, then 



(10.2) 5(0) = ±2 / luiKdC, 

where, in the left hand side, one integrates in the increasing direction on the real axis. 

Proof. Deform the integration contour G{q) so that it go around q just along it. Then, relation (10.2) 
follows directly from (10.1). □ 

This lemma immediately implies 

Corollary 10.1. Let E £ J. Then, 

(1) S{q) is real and non-zero; 

(2) as a functional of the branch k, it takes only two values that are of opposite sign. 

Proof. Inside any spectral gap, the imaginary part of no branch of the Bloch quasi-momentum 
vanishes. Hence, the first statement follows from (10.2). The second one follows from (4.5) listing all 
the branches continuous on the integration contour. □ 

In the sequel, we choose the branch k so that, on J, S{q) be positive. S{g) is called the tunneling 
action. 

10.4. Obtaining estimates (1.13). All the estimates in (1.13) are obtained in the same way. So, we 
prove only the estimate for b-m in the case of n odd. Recall that we work in Vb, a small constant 
neighborhood of a point £"0 G J. 

The coefficient b-m is given by (9.6). The definition of the amplitude of an arc, formula (8.3), implies 
that A{f, Tf, di) is independent of e, continuous in E and does not vanish. So, there are two positive 
constants Ci and C2 such that 

(10.3) Ci<\A{f,Tf,do)\<C2, E£Vo. 

Let us estimate the factor exp (^5(/, T/, do)) for E £ Vq nW. Therefore, we choose the arc 7 = 
7(7, Tf, do) stretched along the real line and going around the branch points (between ^0 ~ 27r and Co; 
the beginning and the end of 7) along infinitesimally small circles. We compute 



(10.4) 



exp ^^S'(/,r/, do)^ ^^"^^("e/ ^^^9*^^^ = exp y-^^ j ImK^dC 



where Q consists of all the connected components of G H M between ^0 — 27r and Co, i-e. of all the 
intervals — 27r, the interval gn+m — 27r and the interval gn-i- 

Using (9.4), one easily checks that, in the right hand side of (10.4), Im Kg < inside each of the 
intervals of integration (which are segments of the arc 7) . 
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Due to the periodicity of W, we can write 
(10.5) 



gee ■ 

where, on each interval of integration, k is any continuous branch of the complex momentum such that 
Im K < 0. By means of Lemma 10.4, we check that, up to the sign, the expression —2 J Im/td^ is equal 



exp ( ^S{f,Tf,do] 



to S{q), the tunneling action. As both are positive, they coincide. Therefore, 
JJ(t(g))~^. This and (10.3) imply the estimate for announced in (1.13). 
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